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ABSTRACT 

White dwarfs are the end product of the lifes of intermediate- and low-mass stars and their evo- 
lution is described as a simple cooling process. Recently, it has been possible to determine with 
an unprecedented precision their luminosity function, that is, the number of stars per unit volume 
and luminosity interval. We show here that the shape of the bright branch of this function is only 
sensitive to the averaged cooling rate of white dwarfs and we propose to use this property to check 
the possible existence of axions, a proposed but not yet detected weakly interacting particle. Our 
results indicate that the inclusion of the emission of axions in the evolutionary models of white dwarfs 
noticeably improves the agreement between the theoretical calculations and the observational white 
dwarf luminosity function. The best fit is obtained for m a cos 2 /3 ~ 5 meV, where m a is the mass of 
the axion and cos 2 (3 is a free parameter. We also show that values larger than 10 meV are clearly 
excluded. The existing theoretical and observational uncertainties do not yet allow the confirmation 
of the existence of axions, but our results clearly show that if their mass is of the order of few meV, 
the white dwarf luminosity function is sensitive enough to detect their existence. 
Subject headings: elementary particles — stars: luminosity function, mass function — white dwarfs 



1. INTRODUCTION 

One solution to the strong CP problem of quantum 
chromodynamics is the Peccei-Quinn symmetry (Pec- 
cei & Quinn 1977a, 1977b). This symmetry is spon- 
taneously broken at an energy scale that gives rise to 
the formatio n of a light ps e udo-scalar par ticle named 
the "axion" (jWeinbergj|197l IWilczeklll978l ). This scale 
of energies is not defined by the theory but it has to 
be well above the electroweak scale to ensure that the 
coupling between axions and matter is weak enough to 
account for the lack of a positive detection up to now. 
The mass of axions and the energy scale are related by 
m a w 0.6(10 7 GeV/fg) eV. A strophysical and cosmologi- 
cal arguments (Raffclt 2007) have been used to constrain 
this mass to the range 10 -2 eV > m a > 10~ 4 eV. For 
this mass range, axions can escape from stars and act as 
a sink of energy. 

White dwarfs are the final evolutionary phase of low- 
and intermediate-mass stars {M < 10 ± 2 M Q ). Since 
they are degenerate objects, they cannot obtain energy 
from thermonuclear reactions and their evolution can 
be described just as a gravothermal process of cool- 
ing. Therefore, if axions exist, the properties of these 
stars would be noticeably perturbed. Furthermore, white 
dwarfs have a relatively simple structure: a degenerate 
core that contains the bulk of the mass and acts as an 
energy reservoir and a partially degenerate envelope that 
controls the energy outflow. The vast majority of white 
dwarfs have masses in the range 0.4 < M/Mq < 1.05 — 
although these figures are still uncertain — and have a 
core made of a mixture of carbon and oxygen. All of them 
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are surrounded by a thin helium layer, with a mass rang- 
ing from 10~ 2 to 10 -4 Mq which, in turn, is surrounded 
by an even thinner layer of hydrogen with a mass be- 
tween 10~ 4 and 10~ 15 M , although about 25% of white 
dwarfs do not have hydrogen atmospheres. White dwarfs 
displaying hydrogen in their spectra are called DA and 
the remaining ones are known as non-DAs. Because of 
the different opacities, DA white dwarfs cool more slowly 
than the non-DA ones. 

The standard theory of white dwarf cooling can be 
summarized as follows ()Isern et~allfl99l . When the lu- 
minosity is large, -Mbd < 8, the evolution is dominated 
by neutrino emission. In this phase the main uncertain- 
ties come from our poor knowledge of the initial condi- 
tions. Fortunately, it has been shown that all the ini- 
tial thermal structures conver ge toward a unique one 
(jD'Antonna fc Mazzitellil 119891 ). For smaller luminosi- 
ties, 8 < Mboi < 12, the main source of energy is of 
gravothermal origin. In this phase, the Coulomb plasma 
coupling parameter is not large and the cooling can be ac- 
curately described. Furthermore, the energy flux through 
the envelope is controlled by a thick nondegenerate or 
partially degenerate layer with an opacity dominated by 
hydrogen, when present, and helium, and it is weakly de- 
pendent on the metal content since metals sink towards 
the base of the envelope by gravitationally induced dif- 
fusion. Below these luminosities, white dwarfs evolve 
into a region of densities and temperatures where the 
plasma crystallizes. When this happens, two additional 
sources of energy appear. The first one is the release 
of latent heat during crystallization. The second one is 
the release of gravitational energy induced by phase sep- 
aration of the different chemical species (Garci'a-Berro 
et al. 1988a, 1988b; Isern et al. 1997, 2000). When the 
bulk of the star is solid the white dwarf enters into the 
Debye cooling phase and the only important source of 
energy comes from the compression of the outer layers. 
These late phases of cooling are not yet well understood 
(|Isern et al.lll998f) . 
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2. THE WHITE DWARF LUMINOSITY FUNCTION 

One way to test the evolutionary properties of white 
dwarfs is using their luminosity function, which is de- 
fined as the number of white dwarfs per unit volume and 
magnitude. The fi rst luminosity fu nction was derived 
four decades ago (|Weidemannl [1968) and since then it 
has been noticeably improved. The most recent deter- 
minations use data from t he Sloan Digital S ky Survey 
(SDSS). The first of these (jHarris et alj |20 06 f) was built 
from a sample of 6000 DA and non-DA white dwarfs 
with accurate photometry and proper motions culled 
from the SDSS Data Release 3 and the USNO-B cat- 
alogue, whereas the second one (jDeGennaro et all 2008) 
was constructed from a sample of 3528 spectroscopically 
identified DA white dwarfs from the SDSS Data Release 
4 (see Fig. [IJ. The monotonic behavior of this function 
clearly proves that the evolution of white dwarfs is just 
a cooling process. The sharp cutoff at low luminosities is 
a consequence of the finite age of the Galactic disk. The 
luminosity function of white dwarfs can be computed as: 
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where Mboi is the bolometric magnitude, M is the mass 
of the parent star — for convenience all white dwarfs are 
labeled with the mass of its progenitor — Tq is the age of 
the population under study, i coo i is the time that a white 
dwarf with a progenitor of mass M takes to cool down to 
a bolometric magnitude MboU ''cool is the characteristic 
cooling time of the white dwarf, t coo i = dt C00 \/dMi> \, 
and t ps is the lifetime of the parent star. In equation 
Q M u is the maximum mass of white dwarf progenitors 
and M\ is the mass of the progenitor that satisfies the 
condition Tq — t coo \ — t ps = 0, i.e., the minimum mass 
of a star able to produce a white dwarf of the required 
luminosity. The remaining quantities are the initial mass 
function, <f> (here we have used Salpeter's law), and the 
star formation rate, ij), which is not known a priori and 
depends on the population under study. 

Since neither the star formation rate nor the total den- 
sity of white dwarfs are known, theoretical luminosity 
functions are normalized to a given observational bin, 
usually the one with the smallest error bar. Figure [1] 
displays the observed luminosity functions and several 
luminosity functions obtained with the same DA cooling 
sequences fS alaris et al.l 12000^ and different star forma- 
tion rates and ages of the Galaxy. The cooling models 
assume a nonhomo geneous distribution of carbon and 
oxygen in the core ([Salaris et alJ Il997| ) , a pure helium 
layer of 10~ 2 M* and on top of it a pure hydrogen layer 
of 10 -4 M*, where M* is the mass of the white dwarf. 
Since a relationship connecting the mass of the white 
dwarf and the mass of its progenitor is also necessary 
we have adopted the one tha t best reproduces th e mass 
distribution of white dwarfs ([Catalan et alJ [20081 ). 

An interesting feature of Figure Q] is that the bright 
part of the white dwarf luminosity function — that with 
bolometric magnitude M^ Q \ < 13 — is almost indepen- 
dent of the assumed star formation rate. This can be 
explained with simple arguments. Since the characteris- 
tic cooling time is not strongly dependent on the mass of 
the white dwarf, equation fl} can be written as 



Fig. 1. — Luminosity functions of white dwarfs. Filled squares 
correspond to the luminosity function of white dwarfs (DA and 
non-DA types) obt ained using the reduced proper motion method 
IIHarris ct al. 2006) and open squares to t hat obtained using spec - 
troscopically identified DA white dwarfs liDeGcnnaro et aUI2008l 'l. 
The solid lines are the theoretical luminosity functions obtained for 
different ages of the Galaxy — from left to right: 10, 11, 12 and 
13 Gyr — and a constant star formation rate. The dotted lines 
are the same for an exponentially decreasing star formation rate 
— ip(t) = exp (—t/r) with r = 0.5, 3 and 5 Gyr — and the same 
age of the Galaxy, 11 Gyr. All the luminosity functions have been 
normalized to the same observational data point. 
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Restricting ourselves to bright white dwarfs namely, 
those for which t coo \ is small — the lower limit of the in- 
tegral is satisfied by low-mass stars and, as a consequence 
of the strong dependence of the main-sequence lifetimes 
with mass, it takes a value that is almost independent 
of the luminosity under consideration. Therefore, if ip 
is a well-behaved function and Tq is large enough, the 
integral is not sensitive to the luminosity, its value is ab- 
sorbed by the normalization procedure, and the shape 
of the luminosity function only depends on the averaged 
physical properties of white dwarfs. It is important to 
mention here that the initial-final mass relationship en- 
ters as a weight into the calculation of this average. Nev- 
erheless, since only those functions able to provide a good 
fit to the mass distribution of white dwarfs are accept- 
able, its influence on the shape of the bright branch of the 
luminosity function is minor. Here we use this property 
of the white dwarf luminosity function, together with the 
recently obtained high-precision observational luminosity 
function, to study the influence of the emission of ax- 
ions, and to check which mass of the axion is compatible 
with observations. The idea of using the cooling times of 
whit e dwarfs to co nstrain the properties of axions is not 
new (|Raffeld[l986[ ). but the crudeness of the theoretical 
models and of the observational data prevented a definite 
conclusion and just a loose upper bound was obtained, 
m a < 30 meV. 
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Fig. 2. — Energy losses for a 0.61 Mq white dwarf as a function 
of the bolometric magnitude. The dashed lines represent the axion 
luminosity for different values of m a cos 2 — from top to bottom: 
m a cos 2 /3 =10, 5, 1, 0.1 and 0.01 meV. The thick solid line is 
the photon luminosity, while the thin solid line shows the neutrino 
luminosity. 

3. AXIONS AND THE WHITE DWARF LUMINOSITY 
FUNCTION 

Axions can couple to photons, electrons and nucle- 
ons with a strength that depends on the specific im- 
plementation of the Peccei-Quinn mechanism . The two 
most common imple mentations are the KVSZ (iKimlll97a 
Shifman et all 1 19801) a nd the DFSZ models (jDine et al l 
19811: IZhitnitskiilll980D . In the first case axions couple to 
hadrons and photons, while in the second they also cou- 
ple to charged leptons. For the temperatures and densi- 
ties of the white dwarfs under consideration, only DFSZ 
axions are relevant and in this case they can be emit- 
ted by Compton, pair annihilation and bremsstrahlung 
processes, but only the last mechanism turns out to 
be important. Figure O shows the energy losses for 
an otherwise typical 0.61 Mq white dwarf as a func- 
tion of the bolometric magnitude. The dashed lines 
represent the axion luminosity for different values of 
m a cos 2 (3. The axio n emission rate (in erg g -1 s -1 ) 
has been computed dNakagaw a et al.l 119871 . Il988f) as 
e a = 1.08 x W 23 a{Z 2 /A)TfF, where F is a function 
of the temperature and the density which takes into ac- 
count the properties of the plasma, and a = 5 2 e /47r is 
related to the axion-clcctron coupling constant g ae = 
2.8 x 10~ n m a cos 2 j3/l eV. Since the core is almost 
isothermal it turns out that L a oc T 4 in the region in 
which axions are the dominant sink of ener gy. The thick 
solid line represents the photon luminosity (jSalaris et al.l 
120001) . For the region of interest L v cx T a , with a ~ 2.6, 
although this value changes as the white dwarf cools 
down. The thin solid line represents the neutrino lu- 
minosity, which scales as L u oc T 8 and is also dominated 
by th e plasma and bremsstrahlung processes (|Itoh et al.l 
19961). Therefore, since the temperature dependence of 
the different energy-loss processes is not the same, the 
luminosity function allows to disentangle the different 
contributions. 
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Fig. 3. — White dwarf luminosity functions for different values 
of the axion mass. The luminosity functions have been computed 
assuming m a cos 2 f3 = (solid line), 5 (dashed line) and 10 (dotted 
line) meV. 

Figure [5] shows that in the region Mbd ~ 10 the axion 
luminosity can be comparable with the photon and neu- 
trino ones, depending on the adopted axion mass. It also 
shows that the region around Mbd ~ 12 provides a solid 
anchor point to normalize the luminosity function be- 
cause there the observational data have reasonably small 
error bars, models are reliable, neutrinos are not relevant 
and axions, if they exist, are not dominant. 

4. RESULTS AND CONCLUSIONS 

Figure [3] displays several luminosity functions obtained 
using different axion masses, adopting a constant star 
formation rate and an age of the Galactic disk of 11 
Gyr. As already mentioned, it is important to real- 
ize that the bright branch of the luminosity function is 
not sensitive to these last assumptions. All the lumi- 
nosity functions have been normalized to the luminosity 
bin at log(L/L ) ~ —3 or, equivalently, Mboi — 12.2. 
The best-fit model — namely that which minimizes the 
X 2 test in the region — 1 > log(L/L ) > —3 (that is, 
7.2 < Aiboi < 12.2), which is the region where both the 
observational data and the theoretical models are reliable 
— is obtained for m a cos 2 /3 « 5 meV and solutions with 
m a cos 2 j3 > 10 meV are clearly excluded. Figure 0] dis- 
plays the behavior of % 2 as a function of the mass of the 
axion in our fiducial case (solid line) and in the case in 
which we use the initial-final mass relationship of Wood 
(1992), which is margina lly compatible with the white 
dwarf mass distribution ([Catalan et al.l 12008). In both 
cases the behavior of the luminosity function is similar 
and gives similar values for the mass of the axion once 
one takes into account the present uncertainties. It is 
also important to notice that the largest contribution to 
the lack of accuracy comes from the brightest bins of the 
luminosity function, which have large error bars. 

This result is completel y compatible with the previ- 
ously existing constraints (|Raffeltll2007h . Furthermore, 
these values are also compatible with the bounds im- 
posed by the drift of the pulsational period of the ZZ 
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4. — Value of x 2 as a function of the mass of t he axion for the 
case i n which the initial-fina l mass relationship of Catal an et al.l 
12008ft (solid line) and that of lWoodl ([1992) (dotted line) are used. 
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Ceti star G117-B15A (llsern et al l Il992t ICor sico et al.l 
[2001 Bischo ffTChn et al Il2008t )~ It is worthwhile to men- 
tion here that axions with m a « 5 meV would change 
the expected period drift of variable DB white dwarfs 
- w hich have values between P ~ 10~ 13 and 10~ 14 
s s _1 (j Corsico fe Althausl [20041 ) — by a factor 2, the ex- 
act value depending on the adopted temperature of the 
stellar core and that the detection of such a drift would 
provide a strong additional argument in favor of the ex- 
istence of axions. 

The results presented here are not a definite proof of 
the existence of axions, since there are still some ob- 
servational and theoretical uncertainties. However, the 
calculations reported here show that the hot branch of 
the white dwarf luminosity function is a powerful tool 



to test the existence of weakly interacting massive par- 
ticles because it is only sensitive to the averaged cooling- 
rate of white dwarfs and not to the details of the star 
formation rate or the initial mass function, as shown in 
§2. Moreover, our results are problably the first evidence 
that the shape of the white dwarf luminosity function 
could be affected by the emission of axions and that this 
change can be measured. If this is indeed the case the 
mass of the axion would be of about 5 meV. In addi- 
tion, we have derived an upper bound to the mass of 
the axions of 10 meV, which is compatible wit h other 
recent determinations (Bischoff-Ki m et all l2008t ). It is 
worth mentioning that this result is relevant for other 
research fields and, in particular, for cosmology. Specifi- 
cally, assuming cos 2 (3 = 1 , the contributio n of axions t o 
dark matter would be of the order of 0.2% (|RaffeltJl2007t ). 
In view of this, we consider it of the largest importance 
to improve the observational determination of the white 
dwarf luminosity function, especially in the region of the 
hottest white dwarfs (see Fig. [3|). Thus, the extension 
of the SDSS is of the maximum interest not only for as- 
tronomers and cosmologists, but also for particle physi- 
cists. However, not only observational efforts are needed, 
since in order to obtain a reliable determination of the 
mass of the axion it is also important to decrease the 
uncertainties in the plasmon neutrino emission rates at 
the relevant temperature range. Furthermore, it would 
be also convenient to intensify the study of the drift of 
the pulsational periods of variable white dwarfs in order 
to obtain additional independent evidence. 
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